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Abstract 

The present series of three consecutive papers develops a general theory for linear and finite solid 
viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the 
general approach is specified in a form convenient for solving problems important for many industries that 
involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of 
viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. 

In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is 
applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some 
non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons 
to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear 
viscoelasticity. The most important is the thermodynamic modeling of “thermo-rheological complexity”, 
i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A 
special structure of interaction matrices is established for different physical mechanisms contributed to the 
normal relaxation modes. This structure seems to be in accord with observations, and creates a simple 
mathematical framework for both continuum and molecular theories of the thermo-rheological complex 
relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining 
both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers 
in the rubbery and glassy regions. 
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1. Introduction 


The general mathematical structure of the linear theory of viscoelasticity has been 
established almost fifty years ago by Gross (1953), with more details elaborated by Bland 
(1960). In accordance with these texts, modeling solid-like linear viscoelasticity should 
include a linear hereditary functional with a convex kernel representing the relaxation 
function, and an out-of-integral linear elastic term that models the equilibrium elastic 
effects. In another approach, Meixner (1953,1954) and Biot(1955) developed the theory 
of linear viscoelasticity on the basis of non-equilibrium thermodynamics, whose general 
theoretical framework was established later by de Groot and Mazur(1962). Numerous 
applications of linear viscoelasticity to characterization and analyses of viscoelastic 
properties of liquid and solid polymers were considered in the text by Ferry (1980). If the 
volume (bulk) deformations are also important in linear region, the linear viscoelastic 
description should generally include two viscoelastic spectra, one for shearing and 
another for bulk viscoelasticity. Leonov (1996) recently discussed the non-equilibrium 
thermodynamic approach to non-isothermal bulk relaxations. It is well known Ferry, 
1980) that the bulk relaxations are usually much faster than the shearing ones. Therefore 
except studying solely bulk relaxations, when both the shear and bulk relaxations proceed 
together, the effects of volume relaxations can be ignored. 

The big obstacle in constructing a uniform thermo-mechanical description of 
polymers is their “thermo-rheological complexity”, well documented during the last 
decade in many papers by Plazek, Ngai, Roland and coworkers (e.g. see: Plazek et al, 
1992; Ngai et al, 1997; and Roland et al, 1997). In these papers, different temperature 
dependences for different parts of relaxation spectra were observed over a very broad 
region of relaxation times or frequencies that includes all the basic types (flow, rubbery 
and glassy) of polymer mechanical behavior. 

In the flow and rubbery regions, the uniform Williams-Landell-Ferry (WLF) 
time-temperature scaling (e.g. see Ferry, 1980) is valid for relaxation parameters 
describing different relaxation modes. This scaling uses the concept of entropy elasticity 
(e.g. see Treloar, 1975), where the equilibrium elastic modulus is proportional to the 
absolute temperature, and the same temperature dependencies apply for modal viscosities 
and elastic module. In these regions, the relaxation spectra for many polymers are well 
described by the Prony series with a relatively small number of relaxation modes. 
However, in the glassy region where the fractal properties of polymers dominate, another 
time-temperature scaling exists. The viscoelastic effects in this region are well described 
by a continuous distribution of relaxation spectrum, such as the stretching exponent 
and/or fractional derivative approaches (e.g. see Schiessel et al, 2000). Employing the 
Prony series in this region, though possible, is impractical since the discrete relaxation 
spectrum, being densely populated, would include numerous discrete relaxation modes. 

It has been demonstrated that a possibility to describe the thermo-rheological 
complexity is by matching both the time-temperature scales when using a conceptual, 
semi-empirical, “mode coupling” model (e.g. see Ngai et al, 1997). This model sometime 
uses more than two matching procedures, when additionally the Rouse relaxation region 
(e.g. see Ferry, 1980) is involved. Some additional arguments have been applied within 
context of this model to cross-linked rubbers (Roland et al, 1 997). 
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Several molecular-based theoretical approaches to linear relaxation phenomena in 
cross-linked rubbers have also been recently developed (Vigilis and Heinrich, 1994; 
Bastide and Boue, 1996; Gotlib and Gurtovenko, 1997). In these theories, a rubber has 
been modeled as “molecular network” filled with a Newtonian liquid, with the polymer 
chains in network being modeled as in the classic statistical theory of rubber elasticity. 
The interactions between the liquid and network and between the segments of the same 
strands were modeled as different types of viscous friction. These approaches are aimed 
at a theoretical derivation of discrete relaxation spectra with only a few physical 
parameters. Whereas Gotlib and Gurtovenko (1997) analyzed the homogeneous polymer 
network with a “classical” elastic contribution, Bastide and Boue (1996) studied the 
effect of topology on the network elasticity, and Vigilis and Heinrich (1994) analyzed the 
effects of the network heterogeneity, such as different types of disorder in the network 
structure and topology. 

Some efforts were also made to develop mathematical procedures for finding the 
relaxation spectra directly from experimental data, which has long being established as 
generally an ill-posed problem (e.g. see review in paper by Simhambhatla and Leonov, 
1993). Several regularization methods with the use of computerized procedures have 
been developed in the literature to overcome this difficulty. Nonetheless, these 
procedures of extracting the Prony parameters still depend highly on initial guesses about 
the numbers of relaxation modes, and usually operate with very high numbers of these. 
Fulcheron et al (1993) and Simhambhatla and Leonov (1993) demonstrated that the Pade- 
Laplace method is a more reliable tool for experimental determination of relaxation 
spectra when represented by a Prony series. It operates with an analytical continuation of 
theTelaxation function into the complejr plane whefe the prdblem““df extraction^ of 
relaxation spectra is well posed and describes the Prony series with considerable less 
numbers of relaxation modes. A computerized Pade-Laplace procedure for extracting 
relaxation modes from a combination of different experiments was developed and 
demonstrated in a paper by Simhambhatla and Leonov (1993) and in the Ph.D. thesis by 
Simhambhatla (1994). In practical applications of the procedure to a viscoelastic solid, 
the true elastic (out-of-integral) component should be found first and then excluded from 
the Pade-Laplace procedure. The same is true for the cross-linked polymers or 
elastomers. No quantitative application of such an approach to solid-like viscoelasticity 
has yet been demonstrated. 

The structure of the present paper is as follows. We first consider the general 
thermodynamic scheme for analyzing the solid-like viscoelastic phenomena, much 
similar to that published earlier in papers by Meixner (1953,1954), Biot (1955) and 
Leonov (1996). Then at the end of the paper we discuss the extension of our scheme to 
incorporate the effects of thermo-rheological complexity. We analyze in this paper only 
incompressible (shearing) solid like relaxations. This is because out of two (bulk and 
shearing) independent, linear relaxation phenomena, the non-isothermal bulk relaxations 
have been recently discussed by Leonov (1996). 

2. State Variables and Free Energy Function 

We hypothesize that the complete set of macroscopic state variables, sufficient to 
describe the linear shearing viscoelasticity are: the temperature T , the actual (total) 
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infinitesimal, traceless Hencky strain tensor h Q , and the set of “hidden” tensor variables 
{=*} (£ = 1>2 depending generally on the space-time variables x,t. Although a 

linear theory does not necessitate knowing the physical sense of the hidden variables h _ k , 

we can guess, and justify later, that they represent elastic strain contributions in the total 
viscoelastic strain originated from k-th physical mechanisms. Due to this guess we 
assume that the tensors h_ k are traceless, symmetric, second-order tensors. We also 

assume that the hidden variables h , describing pure non-equilibrium properties of the 
system, vanish in the equilibrium. 

We now introduce the specific (per mass unit) Helmholtz’ free energy 
F depending on the above set of state variables: 

F = F(T\h Q ,h x ,...,h n ) (trh Q =trh k =0). (1) 

The specific entropy, S , and all other specific potentials, such as internal energy U , can 
then be found using common thermodynamic formulas, including, the Gibbs’ relation: 

dU = dF + TdS, U = F + TS, S = -dF / dT^, etc. (2) 

Near the equilibrium, unloading state, which is far away from possible phase transitions, 
the free energy F is represented as a quadratic, positively definite form: 

pAF.W = V2G,(T)tr£+l/2 f l G v (T)lr(h [ -h ,). (3) 

*J= 1 

corresponding to the minimum of the free energy. Here W is called the elastic potential. 
The positive definiteness of form (3) means that the equilibrium shear modulus G 0 is 
positive and the nxn symmetric energetic matrix, ||G^(7’)|| is positive definite, too. 

In thermodynamic equilibrium under a stress action, 


<7 =dF/dh 
=0 =0 



= G 0 h o , dF/dh k 



(4) 


Here a Q is the equilibrium extra-stress tensor, i.e. stress without isotropic pressure. The 
second formula in (4) is a consequence of our assumption that the hidden variables 
h_ k vanish in equilibrium. Namely because of this, the possible interactive term 

^ M k (T)tr(h Q h k ) is absent in (3). 

k = i ~ 
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3. Clausius-Duham Inequality and Non-Isothermal Relaxation Phenomena 

To analyze the behavior of the system in non-equilibrium, the common expression for the 
entropy production (or source of entropy) P s is employed in the form: 


TP S = -T~'q-Vr+tr(a-e)-dW/dt\ T (>0). (5) 

Here q is the thermal flux, V is the space gradient symbol, e is the (traceless) strain rate 
tensor, and a_ is the extra stress tensor. The inequality in the bracket of (5) indicates that 

according to the second law, the entropy production P s must be positive for all non- 
equilibrium processes and vanish in equilibrium. With this inequality, the right-hand side 
of (5) represents the well-known Clausius-Duham inequality. 

Substituting (3) into (5) yields: 

TP s =-r'q-VT + tr(a-e)-J^tr(g k -dh k /dt) (>0), (6) 


where 

*=i 

Instead of tensors dh /dt,we now formally introduce new kinematical variables, 
the second-order traceless tensors e p , satisfying the “evolution equations”: 

dh k /dt + e p k =e (k = 0,1, 2,. ..,«). (8) 

Evidently e p =0 for k = 0 , since the well-known kinematical relation for h Q is: 

dh o /dt = e . (9) 

For k > 1 , the unknown tensors e p in (8) characterize the irreversible effects, and in fact, 

the evolution equations (8) simply represent the definition of tensors e p . 

Substituting now (8) with account for (9) into (6) yields: 

TP s =-T-'q-Vr + tr[(g-g T )-e] + ^tr(g k -e p ) (>0), Z T =1tZk ' ( 10 ) 

k = 1 k = 0 

Here q_ T is called the thermodynamic (extra) stress tensor, since it is completely 
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determined by the free energy function W . 

Equation (10) distinctly demonstrates three sources of dissipation in the system: 
(i) non-isothermality, (ii) the deviation of the actual stress a; from the thermodynamic 

one , and (iii) non-equilibrium dissipative effects originated in the above n irreversible 

physical mechanisms. The last two terms in (10) are called the mechanical dissipation. 
The sense of the second source of irreversibility is easily demonstrated in the case when 

the irreversible contributions in viscoelasticity, the tensors h _ k , vanish. Then <x r = !Z 0 > 

and the second term in the first equation in (10) correspond to a dissipative mechanism of 
the Kelvin- Voigt type. This mechanism denies “instantaneous elasticity”, aas observed in 
experimental studies. The existence of instantaneous elasticity necessitates & T =a Q . 

Thus for a modeling of the observable viscoelasticity with instantaneous response, (10) 
takes the final form: 


TP s =-T-' r VT+f d tr(Z t -e P i ) ( 20 ), £ = £ r =Z&- 00 

Jt=l *=0 

The entropy production in (10) (as well as in (11)) is represented as a typical 
bilinear form of the type ^X m -Y m . The independent quantities VT and q_ in (1 1) are 

m 

treated as the “thermodynamic forces”, X m , and conjugated to them are the independent 

quantities, q and e p , known as the “thermodynamic fluxes”, Y m . Note that the 

thermodynamic fluxes are represented via time derivatives that change their signs under 
the reverse time transformation, t-±-t. In the linear scheme of non-equilibrium 
thermodynamics, the independent thermodynamic forces and fluxes are related via linear 
phenomenological relations with temperature dependent scalar coefficients. These 
relations include, in particular, the constitutive equations for the relaxation phenomena. 
Although the general case (10) is easy to analyze, we consider below only more 
physically meaningful and simple case described by (11). Then the linear, 
phenomenological, constitutive relations take the form: 


n 


m= 1 

(12) 

q = -l(T)VT. 

(13) 


Such basic principles as simple arguments of the tensor dimensionality (Curie’s 
principle) and Onsager’s principle of symmetry for kinetic coefficients have been 
taken into account in the constitutive relations (12) and (13). Substituting these into (1 1) 
yields: 


77> s =rU|wf+ < s °) 

*J=1 


9L~9L t 


n 



(14) 
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Here the second term in the right-hand side of the first relation in (14) represents the 
contribution in entropy production due to mechanical dissipation. According to the 
second law, the entropy production must be positive in non-equilibrium. Because the two 
terms in (14) are independent, the coefficient of heat conductivity, X , is positive and the 

kinetic matrix ||^ j| is positive definite, too. 

Expressing e p via gr due to (12), and in turn, ^ via /^due to the second relation 
in (7), yields: 

d^/dl + f t A^(Ty3 v <T)h=e,dhJdn S, = a «<- r >b, + i. a ^(T)h (15) 

mj = 1 k,j = 1 

(k = 1,2,..., n) 

Equations (15) represent the general formulation of linear constitutive equations for 
viscoelastic solids with instantaneous elasticity. They also justify our previous guess that 
the hidden variable h is elastic contribution from the k-th non-equilibrium physical 
mechanism. 

4. Normal Relaxation Modes: Thermo-Rheological Simplicity 

In order to demonstrate that the general, linear, solid, viscoelastic constitutive equations 
(15) reduce to the normal (or independent relaxation) modes, we initially employ the 
assumption of the thermo-rheological simplicity: 

G kJ (T) = G 0 ( T)g „ , A h (T) = a „ / 7o (D . (16) 

Here G 0 (T)(~ T) is the shear modulus known in rubber elasticity (e.g. see Treloar, 1975), 
tj 0 (T) is a characteristic viscosity that far away from the glass transition is described by 

the Arrhenius dependence; |g^| and |a^.| are numerical nxn, symmetric and positive 

definite matrices. The application of general scaling formulae (16) to polymeric materials 
has been widely discussed by Ferry (1980). Although the assumption of thermo- 
rheological simplicity has been criticized in cited papers by Plazek, Ngai, Roland et al, it 
still approximately holds for polymers above their glass transition temperature for 
restricted frequency and time intervals. 

Because both the numerical nxn matrices, ||g^.||and ja kj || , are symmetric and 

positive definite, it is possible to introduce the “normal coordinates”, tensors h k (t,x) , as 
follows: 

£=i>*Sj (i7) 
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Here || q u || is a numerical orthogonal nxn matrix that reduces simultaneously both the 

matrices |g # | and |a^|| to the diagonal form with respective principal values g k and a k . 
Under the transformation (17), the constitutive equations (15) take the form: 

dh k / dt+h k /0 k (T) = e = dh o /dt ; g = G 0 (T)h Q +f j G k (T)h k (k = 1,2,...,#*). (18) 

i=l 


Here the parameters G k ( T ) and 9 k (T) are defined as follows: 

G „ CO = G„ (T)g„ ; e t (D = 1 $ 0 (T)stj 0 (T)/G 0 (T). (19) 

a k g k 

Without loss of generality, we now assume that components of the numerical 
matrices ||g^. || and ||a^|| are sorted so that the products of their principal values, g k a k , 

monotonically decrease with increasing number k, i.e. downwards along the main 
diagonal. Therefore, due to (19), 

a kSk — a k +\gk+i’ d k (T) > 0 k+x (T) . (20) 

The choice of such a monotonous series is needed in the discrete case for formulating the 
common presentation of a relaxation spectrum, G = G(0(T)) k , with a possible formal 
continuous generalization. Evidently, the series G k (T) is generally non-monotonous. 

The free energy function (3) and the contribution of mechanical dissipation in the 
entropy production (the second term in the first equation in (14)) can be expressed via the 
normal coordinates as follows: 

W = \l2G 0 {T)H^ + ^/^G k {T)tr(h k ), (3a) 

*=1 

D = TP s \ T ='£[G l ,<X)ie t (T))rr(g) (> 0) . (14a) 

*,7=1 

Consider also the integral form of the above formulas in isothermal case when 
G k (T) and 0 k (T) are constant. Integrating the first equation in (18) and substituting the 
result in the second one, yields: 

£(0= jexp [-{t-T)ld k ]dh Q {r), a(t) = G 0 h f) (t)+ jig(t-T)dh 0 (r)dT . (21) 

—CO —00 


Here y/(t) is the relaxation function defined for the discrete relaxation spectrum 
{G k ,0 k } as the Prony series: 
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( 22 ) 


KO = l]^*exp (-t!9 k ). 

k = 1 

Formulae (3a) and (14a) can also be represented in the form: 


W = l/2G 0 (TMh 2 o (t)\ + l/2j jv'(2t-T l -T 2 Mdh o (r l )-dh o (r 2 )]. (3b) 
D = -1/2 J JV$2r-r, - T 2 )tr[dh o {r x )- dh o (r 2 )]. (14b) 


Here yJ$2t-r x -r 2 ) = dy/(t) I dt\,_+ 2 ,_ Ti _ Ti (<0). Formulae (21), (22), (3b), and (14b) are 

well known in the linear viscoelasticity with a discrete relaxation spectrum (Gross, 1953; 
Bland, 1960). 

In accordance with observations, the temperature variations of the elastic modulus 
G 0 in the non-isothermal case are usually neglected. Then the new time-like, non- 
dimensional variable t * is introduced as follows: 

t 

t*= \drie(T{r)). (23) 

-00 


Using the characteristic modulus G 0 (T) and relaxation time 9 0 (T) as scaling 
parameters, it is possible to represent the constitutive equations (18) in non-dimensional, 
temperature independent form, which demonstrates the well-known time-temperature 
superposition principle. When it is applicable to experimental viscoelastic data, they are 
represented in the form of “master curves”. The dependence G 0 (T) is commonly used in 
the form of that known in the rubber elasticity (Treloar, 1975), whether the solid-like or 
liquid-like viscoelasticity is studied. The dependence Q 0 (T) above the glass transition 
temperature is commonly used in the WLF form, which at higher temperatures is reduced 
to the familiar Arrhenius form. 

5. Normal Relaxation Modes: Thermo-Rheological Complexity 

In this case both thermodynamic matrices, |G Jjt || and \\A ik || can be transformed under 

isothermal conditions to the diagonal form, i.e. the isothermal problem once again can be 
reduced to the analysis of independent, normal relaxation modes. However, the 
factorization (16) is not valid here anymore. It means that there is a dependence of the 
temperature parameters on the number of the discrete relaxation mode. In other words, 
different temperature dependences exist in this case in various parts of relaxation 
spectrum. Luckily, the experimental observations by Plazek, Ngai, Roland and their 
coworkers, cited above, found these temperature dependencies to be relatively simple, 
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with mostly two consequent intervals of the relaxation spectrum having different 
temperature scales. According to the thermodynamic approach demonstrated in this 
paper, it means that there are negligible physical interactions between these almost 
independent parts of the relaxation spectrum. 

This situation found a sound physical interpretation, especially in the case of 
polymer melts with very long and flexible polymer chains and narrowly distributed 
molecular weight (e.g. see Ngai et al, 1997). It has been revealed that for these polymers, 
three different relaxation spectra can be introduced that relate to the three different scales 
of molecular motion. These are: (i) the supra-molecular scale where motions of large 
entangled macromolecules, which being almost independent on the intra-molecular 
chemical chain structure, contribute to the dynamic module in the flow and plateau 
regions; (ii) the intermediate, Rouse’s molecular scale of motion, which is independent of 
the presence of entanglements, and being nearly independent of intra-molecular chemical 
chain structure, contributes to the ascending (after plateau) region in the dynamic 
module; and (iii) the glassy region with a small intra-molecular scale of motion, which 
depends highly on the chemical structure of monomer units in polymer chains. 
Remarkably, the molecular motions within these three scales can be considered as 
practically independent Almost the same physical picture can be found in (usually 
slightly) cross-linked rubbery materials (Roland et al, 1997), except for the first large- 
scale molecular motion region, which is now mostly affected by the cross-links. This is 
the case considered in the present paper. 

In our opinion, the intermediate, Rouse’s temperature scaling should have the 
same WLF temperature dependence of typical relaxation phenomena, as the first large 
scale motions of polymer chains. This scaling has been introduced by Ngai et al (1997) to 
narrow the “disagreement band” and obtain perfectly smooth master curves. Since the 
temperature scaling is in principle different only in two regions, glassy and flow/rubbery, 
we develop below an example using a two-scale approach. However, there is no 
difficulty to extend it for as many temperature dependences as needed. 

In order to describe this complex thermo-reological behavior, we make a 
simplifying assumption that the above two molecular motions are completely 

independent. In our formal scheme it may happen only if both the energetic, G „ = || G ik || , 

and kinetic, A„ = ||4*||> nxn matrices have similar diagonal block structures (see 
Appendix). It means that the suggested structure of these matrices is that of similar 
diagonal block matrices (DB) matrices that represent the above two main molecular 
mechanisms of polymer motions as: 

G„=diag{G rl ,G r2 },A n =diag{A rV A r2 }, G„ ~ A„ (/*, +r 2 = n) (24) 

The corresponding rank of r k component matrices G r> and A n in Eq.(24) are symmetric, 
positive definite, and are assumed to be represented in the partially factorized form: 

G rj =G*(7X , A* =%(7X t (£ = 1,2); G 0 (T)/G : (T) = c = const. (25) 

Here g n and {k = 1,2) are symmetric and positive definite numerical rank of 
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r k matrices. 

We now use the orthogonal DBM, Q„ = diag{ Q r| ,Q r2 }, (/*, +r 2 = ri) that consists 
of orthogonal rank ^numerical matrices Q of the same structure as shown in Eqs.(24) 
and (25), and reduce the DB matrices G„ and A„ to the diagonal form: 


g; =Q„ G„ q: = Oag{G,<m* gn Qw} 

= diag{G,(T)diag{g k } ri } 


X r < = " 


■£ r .„ — diag{G k (T)} k=X2 n 

(26) 


A* = Q„ • A n • = { 7 , (T)Q ri ■ a„ • Q^} 




= ^{ 7 , ( T)diag{a k }„ } 


= diag{A k (T)} k=xl ,...,n 


It is also possible to introduce normal tensor coordinates, h k (t,x ) , as 

£=1q.iu < 27 > 

that reduce CE’s (15) to the form of independent normal modes: 


* . 2 n+l 

dh l ldt+hje l t (T) = e = dhJdt ■ 2 = G.(r)i 0 +EZ G ;( r &; £'•,="•(28) 

/=! k-ri 


2 

S 

i-1 


Here the parameters G[ ( T ) and 0[ ( T ) are defined as follows: 

G[ CO = G, (Og* ; 0' k (T) = , (29) 

a kSk a kSk 

where 

G, =Gj, 1 < £ < r, ; G t -G 2 , r, <k<r 2 ; (30) 

9i=9\, \<k<,r{, 0i=0 2 , r x <k<r 2 \ . 

We assume without loss of generality once again that inequalities (20) hold. 
Therefore, one can construct the discrete relaxation spectrum G = G(0{T)) k , with an easy 
formal generalization for the continuous relaxation spectrum. Also, in the isothermal 
case, the formulae ( 21 ), ( 22 ) and also (3a,b) and (14a,b) hold true, however, with such a 
modification that the two different parts of the relaxation spectrum ( 22 ) have different 
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temperature dependences. We now assume that 
e,(T)»0 1 (r>. (31) 

These inequalities have been observed in all types of polymers as well as in cross-linked 
rubbers. Then inequalities (20) and (31) guarantee the partial scaling for dynamic and 
relaxation experimental data. 

Consider as an example an interpretation of isothermal dynamic experiments, 
where harmonic oscillations with frequency co are applied to a sample of a cross-linked 
rubber in a very wide range of frequencies and temperature. This is the situation where 
the thermo-rheological complexity cannot be ignored. The expressions for the real 
G'(co) and imaginary G"(co) parts of the complex dynamic module are: 


G\c) = G a + G I f,g„ 
*=1 


W £A) 2 , Q (<»02gk a k) 2 . 

1 + (<o6 x g k a k f 2 * 1 + (c o0 2 g k a k f ’ 


(32) 


G"(©) = G,£ g, 


tl>0\gk a k ! Q f’ ® e 2Sk a k 
\ + {a>9 x g k a k f 1 + (o)0 2 g k a k ) 2 ' 


Using the inequalities (20) and (31), one can show that there are two different 
consecutive “temperature independent” scaling for both G' (tv) and G"(co) in their 
corresponding parts of the frequency range: 


G'M,(r)) G"M(r)) . 

G X (T) ’ G X (T ) 

G’M 2 (7Q) G\(D0 2 (T)) . 

G 2 (T) ’ G 2 (T) 


o <cod x «\l(g rX a rX ) 

(33) 

1 ) < <° e 2 « h • 


Similar types of scaling can be employed to other viscoelastic tests, such as creep or 
relaxation after suddenly applied constant strain, etc. It should be noted that applicability 
of the proposed scaling to creep data is based on the well-known fact that for the discrete 
viscoelastic spectra, the retardation spectrum is “imbedded” into the relaxation one, so 

9 k >9 k > 6 k+l (e.g. see Bland, 1960). Here 6 k are points belonging to retardation 
spectrum. It means that in continuous approximations of the viscoelastic spectra, both 
spectral regions coincide. 

As mentioned, the continuous relaxation spectra are commonly employed in the 
description of relaxation phenomena in the glassy state. This means that the above 
thermodynamic scheme, whic uses a discretized matrix approach with densely populated 
interactive matrices for the fast relaxation phenomena, should be changed to a 
continuous, integral approach. Although this change is quite physically transparent, its 
formal, strict development is still desired. As established here a unified scheme can 
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combine both the long time (discrete) and short time (continuous, or densely populated 
descretized) approaches to relaxations of polymers in rubbery and glassy regions. 


6. Conclusion 

The thermodynamic derivation of linear viscoelastic constitutive equations discussed 
above can be seen as a justification for the completely phenomenological theory of linear 
viscoelasticity. In the thermo-rheological simple case, we used the time-temperature 
superposition principle to reduce the constitutive equations to the form of normal (non- 
interactive) relaxation modes. This derivation is quite similar to the well-known 
derivations from the papers by Meixner (1953,1954) and Biot (1955) (see also Leonov, 
1996). The thermodynamic derivation of linear viscoelastic constitutive equations in the 
thermo-rheological complex case assumes different types of independent equilibria and 
non-equilibrium interactions between the various polymeric subsystems. We have 
modeled these independent interactions using a well-known formalism of diagonal block 
matrices. This approach can also be applied to model dielectric relaxations for thermo- 
electrically complex cases. 

The demonstrated theory creates a solid thermodynamic framework for the more 
detailed, molecular models. These present valuable and complementary parts to 
continuum approaches because they can answer the important unresolved questions about 
the size and “density” of population of interactive elements, along with different 
temperature dependencies of the independent interactions. For example, if the density of 
interactive elements of microstructure is high, as usually happens in the glassy region, 
then enormous amounts of interactive elements should be involved in the rapid relaxation 
response. These rapid relaxations can be successfully described by various continuous 
approaches, such as the “stretched exponential” and fractional calculus approaches 
(Schiessel et al, 2000), related to the fractal physics of these relaxation phenomena. 
Although in this paper we demonstrated the thermodynamic approach as applied to 
discretized relaxation phenomena, the approach could be easily extended to the 
continuous case, when operating with integral equations, instead of matrices. 

The thermodynamic framework for linear viscoelasticity in both thermo- 
rheological simple and complex cases can also be used by applying such a direct and 
mathematically well-posed procedure as developed by the Pade-Laplace method applied 
to viscoelastic phenomena by Fulcheron et al (1993) and Simhambhatla and Leonov 
(1993) for finding the effectively discretized relaxation spectrum directly from 
experimental data. This method can also be applied to different parts of the relaxation 
spectrum under various constant temperatures to reveal their possible different 
temperature dependencies. It is of interest to use the Pade-Laplace method for effective 
discretization of continuous relaxation spectra, such as stretched exponential or trhe 
Mittag-Leffler function,, to evaluate the density of population for the various relaxation 
mechanisms. 

In the case of rubber viscoelasticity, it is also of interest to find the change in 
relaxation spectrum and its temperature dependence during the cure reaction when the 
emerging cross-links change the rubber state from flowing (“green”) to completely cured 
cross-linked rubber. Such an approach has been initiated by Mitra (2000) to reveal the 
change in the linear viscoelastic spectrum during the cure reaction, however, without 
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application of the Pade-Laplace technique. 

Appendix 1: Diagonal Block Matrices 


Let M„ be a nxn matrix and{M r } (r < «)be a set of rxr matrices. We define the 
diagonal block (DB) matrix as a combined M„ matrix representing a diagonal, ordered 
sequence of the component matrices M r : 


M„ =diag{M r[ ,M r2 


..M J= 


M 


rl 


0 

0 

0 


0 0 0 " 

M r2 0 0 

0 0 

0 0 M* , 


'Lh= n - 

k 


(Al) 


Obviously, the order in the diagonal sequence (Al) is essential. 

Consider now the case when the component matrices M r are symmetric and 
positive definite. Then it is evident that for any combination in Eq.(Al), the DB matrix 
M„ is symmetric and positive definite. We call the component matrices properly sorted 
if their principle values increase along the diagonal downwards. A DB matrix is called 
monotonous if the combined principal values of its properly sorted component matrices 
monotonically increase or decrease along the diagonal downwards. The structure of a 
DBmatrix is defined by the consequence of the component matrices of various (in 
particular, the same) allowed ranks. Two DB matrices are called similar if their 
component matrices are of the same structure and have the same properties, i.e. if they 
are, e.g., symmetric or positive definite, or orthogonal, etc. 

Consider now a vector X n e 91 „ and a linear operation in9!„, represented by a 
DB matrix, M„: 


L„=M„ 


■x„=xm y , e *,, 


(A2) 


Equation (A2) demonstrates that any nxn DB matrix M n with a given structure 
fractions the total space 91„ into the full set of mutually orthogonal, non-intersecting sub- 
spaces 91,. of lesser dimensionality. Equation (A2) holds for any DB matrix of various 
structures. For example, if the M r in (A2) are orthogonal, i.e. M,. = Q r , the 
transformation (A.2) is orthogonal. If the M r are symmetric r, xr, component matrices 
in a DB matrix, with principal values m r , and Q r are the r t x r, component matrices of 
the orthogonal DB matrix, which reduce the component matrices M r to the diagonal 
form,M* , then 
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£,.» = d ' 0 «i m »>*-l,!,.J. ' ( A3 ) 


The demonstrated above properties of DB matrices show that they behave quite similarly 
to regular diagonal matrices. 
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